Soliton Junctions 
in the Large Magnetic Flux Limit 



Stefano Bolognesi* and Sven Bjarke Gudnason''' 

The Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 



Abstract 

We study the flux tube junctions in the hmit of large magnetic flux. In this hmit 
the flux tube becomes a wall vortex which is a wall of negligible thickness (compared 
to the radius of the tube) compactifled on a cylinder and stabilized by the flux in- 
side. This wall surface can also assume different shapes that correspond to soliton 
junctions. We can have a flux tube that ends on a wall, a flux tube that ends on a 
monopole and more generic conflgurations containing all three of them. In this paper 
we flnd the differential equations that describe the shape of the wall vortex surface for 
these junctions. We will restrict to the cases of cylindrical symmetry. We also solve 
numerically these differential equations for various kinds of junctions. We flnally flnd 
an interesting relation between soliton junctions and dynamical systems. 
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1 Introduction 



In a recent series of works [1, 2, 3] we studied the behavior of the Abrikosov-Nielsen- 
Olesen (ANO) vortex [4] in the large n hmit, where n is the number of quanta of 
magnetic flux carried by the vortex. We have seen that in this hmit the ANO vortex 
becomes essentially a bag, analogous to the bag models of hadrons [6, 5]. We now 
briefly review the essential results. 

The theory under consideration is the Abelian-Higgs model 

^ = -\f,.F'^'' -\{d,- ieA,)q\' - V{\q\) , (1.1) 

where the potential is chosen such that it has a minimum in the Higgs phase \q\ = 
go 7^ 0. The ANO vortex is a string-like soliton that extends in time and one spatial 
dimension. The simplest way to describe it is to introduce cylindrical coordinates 
{z, r, 9) and orientate it in the z direction. The fields can then be put into the 
following form 

q^qoe'^'q{r) , (1.2) 

Ae = -A{r) . 
er 

The problem is now to evaluate the profile functions q{r) and A{r) subjected to the 
boundary conditions g(0) = 0, g(oo) = q^ and ^4(0) = 0, ^4(00) = 1. The claim is 
that, for every Higgs-like potential V, in the large n limit the profiles become 

lim g(r) ^^H(r-i?v) , (1.3) 
rVi?v^ < r < i?v , 



lim A{r) 



r>RY. 



This conjecture has been proved in [3] by numerical computations. We show in Figure 
1 the result obtained for the BPS potential and n = 25, 000. The step function of 
the profile q{r) reveals the presence of a substructure: a domain wall interpolating 
between the Coulomb phase q = and the Higgs phase q = qo. For this reason we 
have named this object wall vortex. 

The wall vortex is essentially a bag, such as the ones used in the context of the bag 
models of hadrons. A domain wall of thickness Aw and tension Tw is wrapped onto 
a cylinder of radius R. The bag model is a good approximation when the thickness 
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0.25 



0.5 



1.5 



Figure 1: This plot is one of the outcomes of the numerical analysis made in [3]. It shows the 
profile functions q{r) (red/solid line) and A{r) (blue/dashed line) for the BPS potential (/3 = 1) and 
winding number n = 25, 000. 



of the wall is small compared to the radius of the vortex. There are three energy 
terms that must be considered. The first one comes from the tension of the wall and 
is proportional to the radius R. The second comes from the energy density of the 
interior of the bag and is proportional to R^. The third energy term comes from the 
magnetic flux and is proportional to 1/R^. When they are summed we obtain the 
tension as function of the radius R: 

T{R) = + rw27ri? + EottR'' . (1.4) 

The radius of the vortex Ry is the one that minimizes this expression. Physically we 
can understand it in this way. When we derive (1.4) with respect to R we obtain the 
various forces that act on the surface of the bag. The first two terms bring a collapse 
force that tends to squeeze the tube. The third term acts instead as a pressure that 
tends to expand the tube. When the two opposite forces are equilibrated we have a 
stable configuration. 

The aim of the present paper is to study a more generic embedding of the wall 
vortex in the three dimensional space. The configuration (1.4) has the maximal num- 
ber of possible symmetries, it is invariant under cylindrical rotation and z translation. 
Now we want to relax the last condition, namely we consider a wall that has only 
cylindrical symmetry. The radius R is no longer a constant but a function R = f{z), 
and the force balance will not be an algebraic equation but a system of differential 
equations. 
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Solitons in the Higgs phase have received great attention in the last years. A lot 
of different models in different limits have been investigated and the jungle of solitons 
and soliton junctions is enormously vast. Here we just mention the works that have 
mostly influenced the present paper: solitons in nonlinear sigma models [15, 16], 
vortices and walls in supersymmetric models with a Fayet-Ilioupoulos term [7, 8, 9], 
the moduli matrix approach [17, 18, 19], nonabelian vortices and their junction with 
nonabelian monopoles [10, 11, 12], the vortex-monopole- vortex junction [13, 14] and 
also the MQCD realization of soliton junctions [23, 24, 25]. We find particularly 
amazing the various relations between the two fundamental solitons in the Higgs 
phase: the vortex and the wall. More recent works are [20, 21, 22]. The large magnetic 
flux limit is a very general approach, it applies to the Abelian- Higgs model (1.1) that 
is the basic building block of all the more sophisticated theories that contain solitons 
in the Higgs phase. Our result can thus be applied also to these theories and maybe 
explain some of the various relations between walls and vortices. Some of these recent 
developments on solitons are summarized in the reviews [26, 27]. 

The paper is divided into two main parts. In Section 2 we derive the system of 
differential equations that governs the vortex junctions. In Section 3 we explicitly 
study the physical solutions of the differential equations. We conclude in 4 with 
comments and possible applications. 



2 The Differential Equations 



To obtain the differential equations for the profile r = f{z) wc proceed in three steps. 
First in Subsection 2.1 we redo the wall vortex analysis interpreting the minimization 
of T{R) (1.4) as a balance of forces acting on the wall. Then in Subsection 2.2 we 
find the mechanical forces (the ones coming from the tension T-yv and the energy 
density Eq) for a generic profile with cylindrical symmetry. Finally in Subsection 2.3 
we obtain the master equation governing the profile f{z). 
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2.1 Wall vortex redone 



The wall vortex radius is obtained by the minimization of the function T{R). The 
derivative of (1.4) divided by the perimeter 27iR gives 

'"V%+.o^0. (2.1) 



e2i?4 R 

This equation can be interpreted as a balance of forces per unit of area acting on the 
surface of the wall. Preceding from the right to left of (2.1) we have: 



• The energy density is a force per unit area directed inwards; 

• The tension of the wall Tw divided by the radius of curvature i? is a force 
per unit area directed inwards. Note that the other radius of curvature of the 
surface is infinite and does not contribute any extra force; 

• The remaining term — must be interpreted as a force due to the magnetic 
field on the boundary of the cylinder and directed outwards. To check the 
consistency, note that since the flux carried by the wall vortex is $s = and 
the magnetic field is i? = the term is exactly equal to the magnetic 
field energy density 



There is a simple way to understand the origin of the magnetic field force We 



B field 




Ry r 



Fi gure 2: The magnetic field B and the magnetic current J at the boundary of the wall vortex. 
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can write the magnetic field as a function of the radius B{r) where B = inside 
the wall vortex, B — Q outside and the transition is encoded in the function B{r). 
Using the Maxwell equation V A B — J we find J(r) — The force per unit 

— * — * — * 

volume acting on density current is F — J AB thus the force per unit of area acting 
on the surface of the wall is 




and is directed outwards. Note that this force is independent on the precise behavior 
of the function B{r). 



2.2 Mechanical Forces 

Now we want to find the mechanical forces in the case of a generic profile r = f{z)- 
We proceed in this way. First we write the mechanical energy as an integral with 
respect to z of the surface element times the tension Tw plus the volume element 
times the energy density 

^mech = Jdz(2 TtTw / V^TJ^ + TTSq f) . (2.3) 

In absence of the magnetic field, the profile f{z) is the solution to the Euler-Lagrange 
equations obtained by varying E^cch- 

1 _L f 2 _ ft' f 

2nT^^J^j^ + 2ne,f =0. (2.4) 

Dividing by the perimeter 27r/ one obtains the forces per unit area 

1 /" 

^ ^,2)1/2 -^"^{i^ +^0 = (2.5) 

£0 is a force per unit area due to the internal energy density and it does not depend 
on the profile. The other two terms can be identified with the tension divided by the 
two radia of curvature of the surface. To verify it we now compute the two radia of 
curvature geometrically. Take a point P = (r, (p, z) = {f{z), 0, z), or in the Cartesian 
coordinates, (x, y, z) = {f{z), 0, z) on the surface. The tangent vector in that point is 

i^( , ,0, , ) • (2.6) 
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Figure 3: Geometry of the wall vortex surface. 

The plane orthogonal to i and passing through P is 

f (X-/) + (Z-^)=0. (2.7) 

The unit orthonormal vectors lying on the plane (2.7) and orthogonal to i are (j, k), 
where 

= ' ~^ 0, J^-— ) , k = iAj = (0,l,0) . (2.8) 



One radius of curvature is in the plane (i, j) and is the easiest one to compute 

1 _ r 

Rl (1 + /'2)2/3 ■ 



(2.9) 



It is directed outwards in case of positive /" or inwards in case of negative /". The 
other radius of curvature is in the plane (j, k) defined by Eq. (2.7). The curvature 
inwards that lies in the plane (j, k) is 

(2.10) 



R2 f^i + r' 

In summary, the force inwards and outwards at the point P = {f{z), z) is respectively 

1 /" 

^i = y^7^^^w; ^^ = -{1 + fnyr, , (2.11) 

in agreement with Eq. (2.5). 
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2.3 Master equations 



Now we finally write the differential equations that define the wall vortex profile 
r = f{z). First of all we have the Maxwell equations V ■ i? = and V A i? = inside 
the profile. We can rewrite them introducing the magnetic scalar potential defined 
hy B = V(/9. In this way the second Maxwell equation is automatically satisfied 
while the first becomes the Laplace equation for the magnetic scalar potential. The 
boundary conditions are that the B field is tangent to the surface of the wall and 
parallel to the vector i defined in Eq. (2.6). In this way the magnetic flux is constant 
along z and is equal to In summary 

A(p = , V<^ 11 i , $B = . (2.12) 

e 

The other equation is the balance of the forces acting on the wall, that is the gener- 
alization of (2.1) in the case of generic cylindrically symmetric surface: 



B'^ 1 _l_ f '2 _ f" f 



2 



where by B'^/2\^^-^-^ we mean the magnetic field force evaluated at the wall surface. 
Equations (2.12) and (2.13) are a system of coupled differential equations (one partial 
and one ordinary) that defines the profile f{z). 



3 The Junctions 

In this section we will discuss the physical solutions to the master equation (the 
system (2.12) and (2.13)). We first make a rescaling in order to simplify the master 
equation. We rescale the lengths of a quantity / Ryl such that the radius of the 
vortex now is 1. Note that this rescaling applies to all the lengths, so / — > Ryf, 
f f and /" ^ i?vV"- We rescale also the magnetic field such that B = 1 in the 
case of the wall vortex. Thus the rescaled master equations are: 

A(^ = 0, Vv^lli, $B = 7r; (3.1) 
-B\^n-^p )\fl~,,fj + {l-p) = 0. (3.2) 
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It depends only on one parameter p defined by 




(3.3) 



1 - p R^jEq 



The parameter p can take on values from to 1. If p = the collapse force is due 
only to the internal energy density (MIT bag), while if p = 1 it is due only to the 
tension of the wall (SLAC bag). 

3.1 "Near vortex" approximation 

To determine the solution we start with an approximation that is valid when the 
junction is near the vortex (soon it will be clear what we mean). The difficult part of 
the master equations is the partial differential equation that determines the magnetic 
field. To overcome this problem we use an approximation that is valid only when the 
profile is near to be a cone. To determine the B field we proceed with the following 
steps. We take a point P = {z,r = f{z)) on the surface and then we consider the 
cone that passes through it and is tangent to the surface (see Figure 4). Then we 



find the magnetic field that is generated by a charge c at the tip of the cone 




z 



Figure 4: The "near vortex" approximation. 




(3.4) 
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where d is the distance from a generic point {z, r) on the red surface to the tip of the 
cone. In order to have a constant magnetic flux as we vary z, the charge c must be a 
function of z. The magnetic flux is 

9b= f B, = 2-K f drrc(!j) ^—^^ . (3.6) 

Performing the integral we finally obtain: 

- - tttf) ■ 

This expression gives us the function c{z). The magnetic field in the point P is what 
we need for the differential equation (3.2). To obtain it we just take (3.4) evaluated 
at r = f{z)- The function c{z) is given by (3.6) where we have also to remember that 
in our rescaled unit where the radius of the vortex is 1 and the magnetic field inside 
is 1, the magnetic fiux is vr. Finally the differential equation (3.2) becomes: 

W4 1 + f '2 - f" f 

+ P .mJ + (l-P)-0- (3-7) 



(1 + n (yiTr - 1)' + 

As we mentioned before this differential equation is only an approximation. The 
procedure we have used to compute the magnetic field is not exact but becomes 
reliable when the surface is very well approximated by a cone. This means that the 
quantity //" must be very small. Near the vortex the second derivative /" is very 
small and thus we can use the ordinary differential equation (3.7) to understand the 
physical properties of the junctions near the vortex. 

Now we take a look at the solutions to the differential equation. From now on 
we use p — 1 that means only tension and zero energy density. In Figure 5 we have 
the two fundamental junctions: (A) is a vortex connected to a domain wall and (B) 
a vortex that ends on a point. Since this point must be a source of magnetic fiux, 
we can call it a monopole. In Figure 6 we have more general junctions: (C) is a 
wall- vortex- wall, (D) is a monopole-vortex-wall, (F) monopole-vortex-monopole. 

The overall picture becomes more clear if we consider the differential equation 
(3.7) as a dynamical system. This means simply that we rewrite the second order 
differential equation in the following way 

^' N ' (3-8) 

9 = ^{f, 9) 
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Figure 5: The orbits. 

where J-'{f,g) is just what remains if we isolate /" in the differential equation (3.7). 
The differential equations (3.8) are a simple example of a dynamical system. We have 
a phase space (/, g) and a flow deflned on it. The time of the flow is in this case the z 
coordinate. The vortex is the flxed point of the flow / = 1 and g = 0. To understand 
the physical behavior around the flxed point we make a Taylor expansion 




The diagonalization leads to 

± V3 (/ - 1))' = ±v^ ± V3 (/ - 1)) + . . . (3.10) 

What we read from this is the following. The vortex is a stationary point of the 
dynamical system. The linear expansion shows that this is a saddle point. The 
solutions of Figure 5 (with the corresponding z —z reflected ones) are the orbits 
that at z ^ —oo {z — > oo) go into the vortex. The orbits around the vortex are the 
two lines (3.10) g = ±\/3 (/ — 1). If we draw the solutions in a phase plot we obtain 
Figure 7. 

3.2 Flux tube, domain wall and monopole {p = 1) 

The flrst junction we want to study is the flux tube that ends on a domain wall. In 
this case we have to take p = 1 so the Higgs phase and the Coulomb phase are both 
true vacua of the theory. 
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Wall-Vortex-Wall Junction 



Monopole-Vortex-Wall Junction 




Figure 6: More general junctions. 

First we want determine the coefficient of the logarithmic deformation of the wall 
when a vortex is attached to it. What we are going to describe is summarized in 
Figure 8. The piece of junction in the figure has two forces that acts on it. The ffist 
is the vortex that pulls it down with a tension Ty- The second is the wall that pulls 
it up. These two forces must be equal: 

Tv = 27rrrw . ^-^ = , (3.11) 

VI + tanr 

where tan^ = dr/dz. At the end we obtain a differential equation for the proffie r{z): 

r\zf - {2'KTy^/TYfrizf + 1 = 0. (3.12) 

This equation can be trusted only at large r. In fact in making the balance of the forces 
we have not considered the magnetic field inside the wall. For large r its contribution 
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Figure 7: The phase diagram (/,/') around the stationary point (1,0). 




^-^ 2nrT\ycos0 

Figure 8: Asymptotic behavior of the logarithmical bending. 

is negligible since B fall offs as which implies that the energy density B'^r'^ falls 
of as Anyway (3.12) gives the correct coefficient of the logarithmic bending at 
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large r: 

T 

^ ~ — — logr + cost . (3.13) 

ZTTi w 

Since — f-Rv, and we are choosing units in which R\ — 1, this means that the 
vortex-wall junction asymptotically goes like 

r oc e^^ . (3.14) 



Now it is time to finally compute the vortcx-wall junction. The strategy we use 
is the following. Eqs. (3.10) and (3.14) give us the asymptotic behavior of the profile 
function at z — oo (near the vortex) and z +oo (near the wall). We construct 
a trial function with the desired asymptotic behaviors and with a certain number of 
free parameters. These parameters will be adjusted by our numerical code in order to 
approximate in the best possible way the "real" vortex-wall solution. The more the 
free parameters are the bigger is the space of profile functions that our trial function 
can span. For any given choice of the parameters we have a certain profile function 
and we can thus solve the Laplace equation (2.12) and the mechanical forces. To 
solve the Laplace equation we use a finite element method routine. Once we have the 
mechanical force and the magnetic force for a given profile we can compute the total 
force that is just the sum of the two. For the "real" vortex- wall junction the two forces 
must be exactly equal but opposite in direction in every point on the profile (2.13). 
For our trial function we thus compute the norm of the total force (the integral of the 
total force squared) and then minimize it. The minimization gives, within the space 
spanned by our trial functions, the best approximation to the real junction. For the 
vortex-wall junction the result of the computation is given in Figure 9. We show both 
the junction and the force diagram. We have used the same strategy to compute the 
vortex-monopole junction. The result is given in Figure 10. We can finally plot the 
global phase diagram in Figure 11. 



3.3 Flux tube, domain wall and monopole (p < 1) 

Now we want to study the same junctions in the case p < 1. We can expand the 
"near vortex" approximation (3.7) around the stationary point (1, 0) and we obtain 




(3.15) 



Vortex-wall junction 






Figure 10: Vortex-monopole junction and its force plot (p = 1). 



The diagonalization leads to 

And so the orbits around the vortex are the two hues g = ±^J^^ {f — 1). 

The first junction we want to study is the vortex-walL In the case p < 1 the energy 
density of the Coulomb vacuum is not zero and this means that a stable domain wall 
does not exist. On the other hand we have a beautiful solution for the vortex- wall 
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Figure 11: Phase diagram for p — I. The dashed Hnes correspond the near vortex approximation 
given in Figure 7. 

junction in the case p < 1 in Figure 12. What happens to this solution when p is 
decreased? The fact is that a vortex-wall junction can exist if the two dimensional 
space orthogonal to the vortex is compactified. If we compactify the two dimensional 
space on a circle then its radius Ruax is given by 

= 7j% + Tw27r/?v + eoTiR^ = + ^ovri?^^ . (3.17) 

The previous equation is simply the balance of forces. The tension of the vortex, 
that is equal to the sum of the magnetic field, wall and bulk energies, must be equal 
to the tension in the Coulomb vacuum that is given only by the magnetic field and 
the energy density. With our units ($_b = 7r,Tw = f,£:o = ^-^) the radius Ruax is 
determined by the following equation 

^ (l + ^) = 7^ + ^vri^Lx • (3.18) 

For the vortex- wall junction the result of the computation is given in Figure 12. 

The last junction is the vortex-monopole junction for the case p < 1. This junction 
is the most difficult to compute and is also the one where we have reached least 
precision. The difficult part comes from the shape of the wall vortex surface around 
the point where it intersects the axial line r = 0. Due to symmetry reasons the 
magnetic field in this point must be zero so the total mechanical force must be zero. 
In the case p = 1 the Coulomb energy density is zero and so also the curvature 
of the wall vortex surface must be zero (see Figure 10). When the Coulomb energy 
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Vortex-wall junction in compact space withp=0.9 




6 B 10 12 14 16 




— magnetic force 

— mechanical force 

— total force 



Figure 12: Vortex-wall junction and its force plot (p — 0.9). 

density Eq is not zero and it exerts a force that tends to retreat the surface. The two 
radia of curvature are equal, directed outwards and with modulus given by 



-KCurv 



that in our units means 



Curv 



2p 
1-p 



(3.19) 



(3.20) 



For the vortex-monopole junction the result is given in Figure 13. We also plot the 



Vortex-monopole junction with 




- magnetic force 

- mechanical force 

- total force 
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Figure 13: Vortex-monopole junction and its force plot (p = 0.2). 
global phase for p = 0.2 in Figure 14. 
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Figure 14: Phase diagram for p = 0.2. 

4 Conclusions and Discussions 

We conclude the paper with a discussion on possible generalizations and physical 
applications of the junctions of large n vortices. 

4.1 Soliton junctions and dynamical systems 

One interesting outcome of this paper is the relation between soliton junctions and a 
dynamical system (see [33] for an introduction to the subject). Considering z as the 
time of the dynamical system, the flux tube is a stationary point of the differential 
equation. The linear expansion around it shows that this is a saddle point where two 
lines of orbits intersect. Now think about the junction with a vortex in the middle, 
for example the wall- vortex- wall or the monopole-vortex-monopole. We can have a 
lot of this kind of junctions since the length of the vortex in the middle can be fixed at 
pleasure. From the dynamical system point of view this is very simple to understand. 
Since the vortex is a saddle point we can have orbits that goes arbitrarily near to the 
vortex and then escape. 

In the analysis of Subsection 3.1 we have a very simple dynamical system. The 
phase space is two dimensional and consist of the profile f{z) and its derivative f'{z). 
Now suppose we are dealing with a non-large n vortex. The phase space will be 
an infinite dimensional space This is the space of the field configurations q{r, z) and 
Ai{r, z) and their derivatives. Even in this infinite dimensional space we can think of 
the vortex as a stationary saddle point of the dynamical system. The fact that the 
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vortex is a saddle point is just a consequence of the z ^ —z reflection symmetry. If 
we make a linearization around a stationary point and we have a subspace H^"^ of 
orbits that go to the vortex solution at z ^ — oo, we also have a mirror symmetric 
subspace of orbits that go to the vortex solution at 2; — > +00. 

The theory of dynamical systems is a very vast and cvolutcs branches of mathe- 
matics. What we have used in this paper is just a tiny bit of it. A lot of interesting 
phenomena arise in the study of the global properties of the dynamical systems. It 
would be interesting to adopt this point of view in the study of soliton junctions in 
more sophisticated theories. 



4.2 Web of flux tubes 



It is possible to relax the condition of cylindrical symmetry and look for generic 
stationary configurations of the wall vortex surface. Consider first the case of the 
junction between three flux tubes. A vortex that carries n units of flux can be split 
into two vortices with fluxes rii and n2 where n = ni+n2. As long as p is different from 
zero, the vortices are of type I and this means that the tensions satisfy the inequality 
Ti + T2 > T. A non-trivial junction between the three vortices is thus possible and 
the angles are uniquely determined by the tensions. If p = the angle between the 
two smaller vortices is zero and so the junction is trivial. It is nevertheless possible 
to obtain a four-vortex junction. With these basic junctions we can construct a web 
of flux tubes in three dimensions. 

The master equations (2.12) and (2.13) can be generalized relaxing the condition 
of cylindrical symmetry. We have to take a map from a generic punctured oriented 
two dimensional surface to the three dimensional space. The punctures on the surface 
are mapped to the flux tubes at inflnity and every handle correspond to an additional 
flux tube in the web. The master equation is now a system of two partial differential 
equations. The flrst one is just the Laplace equation (2.12) for magnetic scalar po- 
tential inside the tubes (with the correct boundary conditions) and the other one is 
the balance of forces 



£2 



wall 



+ ^ + ^0 
-Kl -ft2 



0, 



(4.1) 



where now Ri and R2 are the two radia of curvature of the surface. 
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4.3 Confining strings 



A possible physical application of large n vortices is in the context of confining strings 
in large N gauge theories (see [28] and [30] for reviews). We briefly review the hy- 
pothesis pointed out in [2] which is the following. Consider an SU{N) pure gauge 
theory and denote by /c-string the string that confines two heavy probes quark and 
anti-quark in the A;-index antisymmetric representation. We then consider the satu- 
ration limit that is the limit where we send both k and to infinity keeping fixed the 
ratio X = k/N. The large N limit is as usual accompanied by the 't Hooft rescaling 
of the coupling constant = g^N. A useful physical quantity is the ratio of string 
tensions divided by N 

TZ{x, N) is the a quantity with a smooth saturation limit. Now comes an assumption. 
Inside the su{N) Lie algebra there is a particular generator, that up to gauge invari- 
ance is unique, that exponentiated passes trough all the elements of the center of the 
gauge group. We assume that the /c-strings are a sort of dual vortices of this U{1). 
This is of course a non-proved assumption but the nice thing, as we are going to see, 
is that it has a very clear signal that can be tested with lattice computations [29]. 
The fact is that the string tension must now satisfy two constraints. The first is that 
in the free string limit {k fixed while N goes to infinity) the tension must be linear in 
k plus subleading corrections. On the other hand, based on the assumption we just 
made, the tension for the dual C/(l) vortex must also be linear when k is large. The 
only reasonable way to combine these two limits is that TZ{x, N) in the saturation 
limit is the triangular function plus subleading corrections. 

n{x, N) = min {x, 1 - x)+0{l/N) . (4.3) 

If this turns out to be correct, the vortex-monopole junction studied in the present 
paper will be a good description of the /c-string quark junction and will maybe be 
visible in lattice computations with a large number of colors. 
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